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Abstract: We apply the previously proposed scheme of approximately self-consistent 
hard-thermal-loop resummations in the entropy of high-temperature QCD to = 4 su- 
persymmetric Yang-Mills (SYM) theories and compare with a (uniquely determined) -R[4,4] 
• Pade approximant that interpolates accurately between the known perturbative result and 

the next-to-leading order strong-coupling result obtained from AdS/CFT correspondence. 
^ , We find good agreement up to couplings where the entropy has dropped to about 85% 

^ I of the Stefan-Boltzmann value. This is precisely the regime which in purely gluonic QCD 

corresponds to temperatures above 2.5 times the deconfinement temperature and for which 
this method of hard-thermal- loop resummation has given similar good agreement with lat- 
tice QCD results. This suggests that in this regime the entropy of both QCD and AA = 4 
SYM is dominated by effectively weakly coupled hard-thermal-loop quasiparticle degrees 
of freedom. In AA = 4 SYM, strong-coupling contributions to the thermodynamic potential 
take over when the entropy drops below 85% of the Stefan-Boltzmann value. 
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1. Introduction 

Much theoretical effort is presently devoted to the study of QCD in a regime where it 
appears strongly coupled. This effort is mainly triggered by the beautiful data obtained 
at the Relativistic Heavy Ion Collider, which reveal that matter produced in high energy 
nucleus-nucleus collisions is strongly interacting, and may, to a good approximation, be 
described as an ideal fluid with very low ratio of shear viscosity over entropy density rj/S 
(see, e.g., Q and references therein). 

At the same time, new techniques have become available that allow the treatment of 
strongly coupled gauge theories. In particular, J\f = A supersymmetric Yang-Mills theories 
(SYM) in the limit of large number of colors A'^ and strong 't Hooft coupling A = g^N is 
now often taken as a model for hot QCD. Through the AdS/CFT conjecture |2|, ^ one 
can make predictions for the otherwise inaccessible strong-coupling regime, in particular 
real-time quantities (such as transport coefficients) |Q, |5|, §|, where lattice is fraught 
with large uncertainties While r]/S ~ (A^logA)^^ » 1 at weak coupling |jlO|, |ll|, the 



AdS / CFT correspondence gives T]/S = l/47r + 0(A-3/2) at large 't Hooft coupling H, |l|] 



which has been argued to be a universal result in strongly coupled gauge theories with 



gravity duals [T^, |l6| 



In [17 1 the thermodynamic entropy of AA = 4 SYM theory at large 't Hooft coupling 



and large N has been determined from the AdS/CFT correspondence as^ 

5/5„.|(l + l^A-3/H...). (1.1) 

The fact that, in lattice QCD, the value of S/Sq for temperatures above 3Tc appears to 
be comparable to this strong-coupling result for = 4 SYM theory has led to the further 



^In as well as Refs. |l|, H, ^ this result appears with (2A)"^/^ in place of X'^'^ , which 
corresponds to a definition of the coupling that differs from the (universally adopted) standard usage of 
g at weak coupling. The latter requires A = g^N = L'^ /a''^ — '^ngstm-n^N , where L is the curvature scale 
of yldSs, a' the string tension, and ^string the string coupling constant (see also ]2l|]). The comparison of 
weak and strong coupling results in Refs. |l^, (with the exception of is therefore incorrect; 

Ref. does have compatible weak and strong coupling results (apart from mixing A with A^), but quotes 
a then incompatible relation between A and pstring- 
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suggestion that also for such temperatures QCD is still in a strong-coupling regime. While 
in this range of temperature the violation of scale invariance in QCD is negligible (in the 
sense that (e — 3P)/e is fairly small), and therefore does not spoil the direct comparison 
with (conformally invariant) AA = 4 SYM theory, we shall argue that such a comparison 
rather suggests a description in terms of effectively weakly coupled degrees of freedom for 
both theories. 

That such a description can be successful should not be taken as implying that the 
system as a whole is literally weakly coupled (typically A > 1), nor that its properties can 
be obtained through a straightforward expansion in powers of A. In fact, the weak-coupling 
expansion for the thermodynamical potential of = 4 SYM theory is presently known to 
order A^/^ ||, |2|, |o|, 



5/5o = 1-Aa + ^A3/^ + ... (1.2) 

but the extremely poor convergence of this expansions seems to limit its usefulness to very 
small values of the coupling, where the deviation from the ideal-gas result is only of the 
order of a few percent (see Fig. ^ below). In particular, this precludes any comparison 
with the corresponding results at strong coupling. This situation is indeed very much 
like in real QCD, where the perturbative series for the thermodynamic potentials of hot 
QCD 25, |2^, ^ shows convergence only for temperatures more than 10^ times the 
deconfinement temperature where the coupling is so small that deviations from ideal-gas 
thermodynamics are minute. 

However, the same impasse appears also in such simple theories as unbroken massless 
O(A^) (p^ theories, where in the limit N ^ oo the entropy is given by an interaction- free 
expression for quasi-particles with thermal masses determined by a one-loop gap equa- 



tion 1 28 1, while the corresponding perturbative expansion is ill behaved except for very 
small values of the (scalar) 't Hooft coupling. As this example shows, the failure of finite- 
temperature perturbation theory in nonabelian gauge theories is not necessarily due to 
specifically nonabelian nonperturbative effects, but largely due to an incomplete resum- 
mation of screening effects (at least as far as thermodynamical quantities dominated by 
hard degrees of freedom, such as the entropy, are concerned).-^ For scalar theories, a corre- 



sponding reorganization of perturbation theory has been proposed in Refs. |30, 31| using 
simple effective mass terms, and generalized to gauge theories in Refs. |3^ by replac- 
ing the simple mass term by the gauge-invariant hard-thermal-loop (HTL) effective action 



1 34, |3g, ISgl- However, this method generally suffers from artificial ultraviolet divergences 
and moreover does not take into account that the HTL effective action is applicable at 
hard momentum scales only for small virtuality. 



^Of course, nonabelian gauge theories do have specific nonperturbative phenomena such as confinement 
in the chromomagnetostatic sector, characterized by the scale g^T. Depending on the quantity under con- 
sideration, the latter can show up in rather low orders of perturbation theory , but in the thermodynamic 
quantities pressure or entropy, they appear only at 4-loop order, parametrically suppressed by 6 powers 
of the coupling g. We are primarily interested in locating the regime where a weak-coupling expansion 
makes sense (sufficiently above the deconfinement transition in the case of QCD) and where we expect such 
contributions of order g^' to be comparatively small. 
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In Refs. 1 37, three of us have proposed a different scheme to resum HTL effects 

which is manifestly ultraviolet finite and applies the (non-local) HTL corrections only in 
kinematical regimes where these are applicable at weak coupling. When applied to QCD, 
we have found reasonable agreement with available lattice data (see Fig. |l]) for tempera- 
tures down to about 3 times the deconfinement temperature^, suggesting that there the 
entropy of the quark-gluon plasma can be understood in terms of weakly interacting HTL 
quasiparticles, the main effect of the interactions being incorporated into a renormalization 
of the quasiparticle properties. 

In the following, we shall recapitulate our HTL resummation scheme for the entropy 
and extend our results to SYM theories, where they can be compared with the strong- 
coupling result of Eq. (|l]5) . In Ref . , we have recently tested our approach successfully 
in the exactly solvable case of large- A^j gauge theories, which is however essentially Abelian. 
The results from the AdS/CFT correspondence deal with the other extreme of infinite 
number of colors. 



2. Quasiparticle entropy and HTL resummation 

Expressed in terms of (self-consistently) dressed propagators {D and S for bosonic and 
fermionic ones) and dressed self-energies (H and S, respectively), the entropy admits the 
following (exact) representation |^, |39| : 



-2tr 



tr J (27r)4 log - ImURe D 

d^k df{uj 



(27r)4 dT 



Im log5~' - ImSRe5 



+ S' (2.1) 



(n and / are Bose-Einstein and Fermi-Dirac distributions, respectively), where S' is a 3-loop 
order quantity that, loosely speaking, describes residual interactions of the quasiparticles: 
the bulk of the interactions have been incorporated in the spectral data determining the 
properties of these quasiparticles through the dressed propagators. 

In practice, some approximations are needed in order to evaluate the quantities which 
appear in Eq. (2.1). Our approximation scheme is based on a skeleton expansion truncated 



at two loop order, whereby S' = 0. Furthermore, the exact expressions for the self-energies 
that would result from solving self-consistent gap equations are approximated by their hard- 
thermal loop expressions, in a way that is described below. This further approximation 
has the virtue of leading to gauge invariant results, which would otherwise not be the case 
in such an approach. 

In perturbation theory, HTL self energies provide the relevant leading-order corrections 
for soft momenta ~ gT as well as for hard momenta, provided the latter have small 
virtuality. One of the merits of the entropy expression ( p.l[ ) in comparison with the usual 
loop-wise expansion of the pressure (whose first derivative with respect to temperature 
also yields the entropy) is that the contributions of order and g^ are generated by either 



^Similar agreement can be obtained by suitable renormalization-scale optimizations of the three-loop 
result from dimensional reduction 
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soft momenta or hard momenta with smah virtuahty, for which the HTL approximation is 
meaningful. 

Indeed, to leading order in g^, all interaction effects are encoded in the (asymptotic) 
thermal masses of hard excitations and summarized by the simple formula |3^ 

b f 

where the sum is over all bosonic and fermionic degrees of freedom. The asymptotic thermal 
masses can be read from the HTL self energies, for which they provide the scale and which 
read explicitly 

for the spatially transverse and longitudinal gauge bosons and 

2 

^oo(f) =F 1 + k 
1 — loe 



for the parts of the fermionic self energy with chirality equal/opposite to helicity, while for 
scalars there is no momentum dependence, lis = "^^(s)- 

Specializing to unbroken SYM theories (without matter in the fundamental represen- 
tation), for the possible choices of we have a number of adjoint scalar (n^) and fermionic 
(uf) degrees of freedom in addition to the two polarizations of the gauge bosons as given 
by Table I (AA=0 corresponding to non-supersymmetric pure-glue QCD). The asymptotic 
thermal masses of the gauge bosons, adjoint scalars, and adjoint fermions turn out to be 
given by 

2 2 2 + ns + nf/2 ^ 

^ooU) = -^^^ ' (2-7) 

which coincide when supersymmetry is realized (whereas only mao{g) is present anyway 
when AA=0). 
Writing 

A A^/^ 

5/50 = 1 + 02^+03^ + ..., (2.8) 

the coefficient 02 is determined by these results for the asymptotic thermal masses through 
( |2.2| ) as tabulated in Table ^ for the possible values of N . 

At low momenta p < A^/'^T, the degeneracy of the thermal masses of the various 
excitations is broken, but to leading order they are determined by the structure of the 
HTL effective action in a universal fashion. This gives different momentum dependent HTL 
masses for gauge bosons and fermions, while the HTL mass of scalar particles is momentum 
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Table 1: Number of adjoint scalar (n^) and ferniionic (n/) degrees of freedom for different A/", and 
the resulting common asymptotic thermal mass for all excitations. Also given are the effective num- 
bers of freedom as measured by = 5o/(27r^T^/45); the coefficients 02 and 03 in the perturbative 
expression for S/Sq] and the values A = A* where = \S2\- 



independent. For gauge bosons, the static Debye mass differs from the asymptotic thermal 
mass according to m?j^ = 2m?^^^y while the plasma frequency in the gauge boson propagator 
is given by oo^i = I^Ti^j^^^ and the fermionic plasma frequency by Mj = 

Perturbation theory at finite temperature is made ill-behaved by the appearance of 
so-called plasmon terms oc X^^'^, reflecting the existence of collective phenomena 



-hnsm^(^)j- = — m^, (2. 



where Ng = — 1. For the different Af this yields the coefficients 03 in ( |2.8D as tabulated 
in Table |l|. The plasmon term renders S/Sq a nonmonotonic function of A, which after 
going through a minimum attains again the Stefan-Boltzmann value at the uncomfortably 
low values A = A* as given in Table ^ the minimum occurring at A = |A=k. Clearly, strict 
perturbation theory breaks down for such couplings. 

However, as the example of 0(A^ — > 00) A(/)^ theory shows [44|, such a breakdown of 



perturbation theory even occurs in an essentially free theory for quasiparticles (5' = 
in this theory) when the coupling constant dependence of the entropy expression ( p.l| ) 
through thermal masses is expanded out. In Refs. [37, 5^, |3^ we have therefore proposed 



to evaluate the thermal self energies in HTL theory to leading and next-to-leading order 
and to keep the entropy expression in Eq. ( p.l| ) unexpanded in the coupling constant. 

Approximating thus 11 and S in by the leading-order HTL expressions ( |2.3D ~( [2^ ) 
turns out to take into account all of ^2. However, only a quarter of the plasmon term ^3 
is due to the entropy of the soft excitations in the HTL approximation, the remaining 
three quarters arise from next-to-leading order corrections of the asymptotic masses of the 
various excitations. These corrections, which are due to the coupling of the hard degrees 
of freedom with the soft collective excitations, are momentum dependent even for hard 
momenta A; ~ T, but to order A^/^ only their average, weighted according to 

- , f dkkn'(k)Re6U(Lj = k) 

is relevant (with a similar expression for the fermionic case). 

For the minimally supersymmetric case {M =1), the solution to the asymptotic ther- 
mal masses can be taken over from Ref. Eq. (13), by replacing Casimir factors of the 
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fundamental representation by those of the adjoint one. This shows that the degeneracy 
of the asymptotic thermal masses in super symmetric theories carries over to + Sm'^. 
Generalizing to > 1 and 7^ 0, we find 

6ml = -ATmoo^^^l^, (2.11) 

which is valid for all the cases covered by Table I. 

Interpreted strictly perturbatively, this leads to tachyonic masses for A > A,,,. However, 
in scalar 0(A^ 00) models, where the same difficulty occurs, approximating the true gap 
equation by one that is quadratic in m gives surprisingly accurate results even for very 
large coupling. This just corresponds to putting 

2 2 I \^ 2V2 + ns 
"^oo = "IooIhtl - ATmoo ^ , (2.12) 



where the notation ttt-^jIhtl is now used for the leading-order quantities in Eqs. ( |2.6D ~(|27 
and moo is defined as the solution to the above equation.^ 



In Ref. 1 39 1 we have proposed to employ such quadratic gap equations for correcting 



the mass of hard excitations by using ( ^.12 ) as the prefactor in the HTL self energies 



( |2.3| ) and (^^), but only when the momentum is hard as specified by a dividing scale 
A = ^27rTm£)CA- This is to take into account that soft excitations are known to have 
rather different corrections in gauge theories (e.g. the Debye mass squared at next-to- 
leading order receives positive and moreover logarithmically enhanced corrections to order 
A^/^ [^], which equally applies for the adjoint scalars^). The introduction of A means 
that for soft excitations the leading-order HTL masses are kept untouched^. As part of the 
theoretical error estimate, ca is being varied in the range ^ . . . 2. In QCD, an even greater 
uncertainty comes however from the dependence of the coupling on the renormalization 
scale fius whose central value we choose as 2ttT, which is close to the optimal value 



considered in Ref. |45] 



3. Results and discussion 

For comparison with our new results for M = 4 SYM which we display further below, we 
start by exhibiting the results of our corresponding previous calculation in pure glue QCD. 
In Fig. H the result obtained in Ref. is reproduced, where the band marked "NLA" 
corresponds to the HTL entropy evaluated to next-to-leading order according to the above 
procedure with ca = ^ • . • 2 and flus = 2ttT. In the larger band bounded by dashed lines, 

^Since this procedure may look rather ad-hoc at this point, let us emphasize that, in the case of the 
scalar theory at least, it corresponds to a meaningful approximation to the actual gap equation required by 
the condition of self-consistency |3^] . 

^The term proportional to log A in the NLO correction to the adjoint scalars' screening mass is in fact 
identical to that of the Debye mass, which in the static case is likewise carried by an adjoint scalar, namely 

^Repeating the procedure of Ref. [^, the separation scale is only introduced for bosonic modes which 
are solely responsible for the plasmon effect. 



-6- 



S/So Pure glue QCD 

7 6 5 MfiMs='27TT 



. 95 



. 9 



0.85 
0.8 
0.75 



-• 1 2 3 T/T, ^ 

Figure 1: The HTL resummed entropy density of pure glue QCD, divided by its Stefan-Boltzmann 
value as calculated in a next-to-leading approximation (NLA) according to Ref . . The full lines 
are evaluated with renormalization scale /2ms = 27rT and different ca (see text). The upper and 
lower dashed lines correspond to renormalization scales fiMS — 47rT and ttT, respectively. The dark 
band represents the continuum-extrapolated lattice results of Ref. pq| . 




the renormalization scale is also varied, fiMS = . . . 4ttT. (The lines marked ca = 
correspond to also rescaling the HTL masses of soft excitations by the same factor as the 
asymptotic thermal mass, which, as we argued above, does not do justice to the former.) 
These results are compared with the continuum-extrapolated lattice results of Ref. |Q 
where the width of the dark band is meant as a rough indicator of the lattice errors (the 
more recent results of Ref. [^] would be centered about the upper limit of this band). For 
temperatures larger than about 3 times the phase transition temperature one finds a 
remarkably good agreement between our NLA approximation and the lattice results, which 
is perhaps even better than one would have the right to expect (given the relatively large 
values of the corresponding coupling 'constant', ffqco ~ -^Z^' with A|^^g=27rr as indicated 
on the top of Fig. |l|). 

The strong-coupling results provided by AdS/CFT correspondence now give us an 
opportunity to test our procedure in the case of = 4 SYM theory. While no lattice 
results are available, the weak-coupling and the strong-coupling results to order A'^''^ and 
\-3/2^ respectively, together should give an idea of the true result. To make this a bit 
more quantitative, we shall construct a Fade approximant^ which reproduces both, the 
weak-coupling and the strong-coupling results, to the order they are known by a rational 
function of A^^^ or, equivalently A~^/^. The fact that at both zero and infinite coupling a 
finite result is obtained requires that this rational function has polynomials of equal degree 
in numerator and denominator. It turns out that there is just sufficient information to fix 

^Pade approximations have already been considered in Ref. [^^, however only for the weak-coupling 
part. The two variants considered in Ref. turn out to coincide numerically rather closely with the NLA 
results for ca = 1 and 1/2 (see Fig. for these latter results). 
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S/So A/" = 4 super- Yang-Mills 




Figure 2: Weak and strong coupling results for the entropy density oiJ\f — 4 SYM theory together 
with the NLA results obtained in analogy to QCD (cp. Fig. Q), but as a function of A, which here is 
a free parameter. The dashed and full heavy gray lines represent the Fade approximants R[i^i] and 
_R[4 4] which interpolate between weak and strong coupling results to leading and next-to-leading 
orders, respectively. 



uniquely all the coefficients in the ^[4,4] approximant 

_ 1 + aX^/^ + /3A + 7A^/^ + 6X^ 

The fact that the weak coupling result has no term of order A^/^, while the strong- 
coupling result approaches 3/4 at large couplings, where it has no terms of order A^^/^ or 



A ^, completely constrains the coefficients in the denominator of in terms of those in 
the numerator: 

- 4 _ 4 - 4 

The remaining four independent coefficients are then uniquely fixed by the results quoted 
above for the weak-coupling expansion to order A^^^ and, respectively, the strong-coupling 
expansion to order A~^/^, together with the expectation that the subsequent term in the 
latter should be suppressed by at least a factor of A~^. One thus finds 

_ 2(9 + 3^/2 + 7^3) 9 
9^ ' ^"27r2' 

Gratifyingly, the resulting rational function has no pole at positive values of A and it 
interpolates monotonically between zero and infinite coupling, as commonly expected [^. In 
Fig. |2|, the Fade approximant ( ^J] ) to the weak- and strong-coupling results is represented 
by the heavy gray line marked "Fade". In order to get an idea of the convergence of 
successive Fade approximants, the simpler Fade approximant R[i^i] which only reproduces 
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the weak-coupling result to order and the constant limit at infinite coupling is also 
shown (dashed gray line). This comparison suggests that the true result should be not too 
far from the -R[4,4] result, which we thus adopt for assessing the HTL-resummed (NLA) 
entropy discussed above. Evaluating the latter precisely in the same way as we did in QCD 
(but with less uncertainty, since because of the finiteness of the AA = 4 SYM theory there 
is no ambiguity from renormalization-scale dependence), we obtain a result that for^ ca 
between 1/2 and 2 agrees remarkably well with the Pade results up to a coupling A ~ 3, 
whereas the ordinary weak-coupling result to order A^/^ fails for all A > 0.5. 

After this comparison within AA = 4 SYM, we may also try to compare the situation 
between the latter theory and QCD. At the top of Fig. ^ we have indicated the values of 
A = 3g^ for (purely gluonic) QCD, when renormalized with two-loop running coupling at 
the scale p-MS = 27rT and for Tc = 1.14Aqcd- In Refs. |48, ^, 11] (see also Ref. [^]), 



where recently other weak and strong coupling results of AA = 4 SYM theory have been 
compared, it has been argued that values of A in QCD (Aqcd) and in AA = 4 SYM (Asym) 
should be related such that comparable thermal effects like screening masses arise. Now, 
(leading-order perturbative) thermal masses in AA = 4 SYM theory are 6 times as large 
as in purely gluonic QCD, just reflecting the much higher number of degrees of freedom 
for a given color number A^. So, by itself, this comparison would suggest a factor of 6 of 
difference between Aqcd and Asym- However, instead of the asymptotic thermal masses, 
one could also compare the normalization of So, the ratio S2/S0 or 53/52 as done in Table 
I, so that the reduction of Asym over Aqcd could be argued to be anywhere between 1/6 
and 1/1.6. 

In fact, the lattice results for the entropy of (pure glue) QCD and the Fade estimate 
of the entropy of SYM present yet another possibility for such a comparison, and to that 
aim we shall focus our attention at larger temperatures in QCD, where S/Sq flattens out. 
(Close to the phase transition the entropy of pure glue QCD becomes very small, whereas 
5/5o in AA = 4 SYM approaches 3/4 at very large coupling.) In QCD, S/Sq flattens out for 
T > 3Tc, corresponding to Aqcd < 5.3. The pure glue 4D lattice results, which are available 
up to about 5Tc, have S/Sq interpolating between about 0.87 and 0.9 (a few percent larger 
values are obtained in the lattice simulations of Ref. [^^)^. Taking the Fade approximant 
( |3.1| ) for S/Sq as plausible estimate for the true nonperturbative result in AA = 4 SYM, a 
comparison of the two theories then leads one to map the regime T > 3Tc of pure glue QCD 
to couplings Asym < 2.5. At such values neither the (strict) weak-coupling nor the strong- 
coupling expansion is giving reasonable approximations, however our 'NLA' resummation 



Just as in the case of QCD we find that it is indeed crucial that the correction (2.f 1) is apphed only 



for hard excitations, as we have done by introducing a separation of hard and soft physics by the scale 
A = y/2nTmDCA and ca ~ 1, which we take as validation of our procedure proposed in the case of QCD. 
Note that for this separation to make sense it is essential that mo ^ 2'jtT, i.e., A ^ 2n^ ~ 20, which is 
indeed fulfilled in the range where our NLA result is plotted for A/" = 4 SYM. 

®The lattice results for QCD that are frequently referred to are usually results with quarks which have 
not been continuum extrapolated and so typically give somewhat smaller numbers. Reliable continuum 
extrapolations have so far been done only for purely gluonic QCD. According to Ref. the continuum 
extrapolation of the lattice results with quarks can be expected to increase the measured values to bring 
them close to the pure glue ratio. 
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of perturbation theory is remarkably close to the Pade results. (Looking only at the strong- 
coupling result to order A~^/^ one would be led instead to couplings A ~ 5.5 and thus 
closer to a direct identification of Aqcd and Asym- However, for such a value of A, the 
NLO result in the strong-coupling expansion deviates already significantly from the LO 
corresponding result of 3/4, thus suggesting that one should expect large corrections from 
higher orders, as indeed suggested by the Pade results.) 

If = 4 SYM theory really provides some guidance here, this would indicate that, 
above 3Tc, the QCD coupling is not large enough to drive the system into a regime that 
resembles that of the strong-coupling limit of SYM theories. This value of the coupling 
is admittedly too large to validate naive perturbative calculations, but not so large as to 
prevent the simple reorganizations of the perturbative expansion that we have presented 
here, and which seems to imply that the entropy of both theories for the corresponding 
couplings can be accounted for in terms of effectively weakly coupled (HTL) quasiparticle 
degrees of freedom. 

Evidently, there is important additional nonperturbative physics at couplings A > 4 in 
SYM as there is in QCD for T < 2.5Tc, though in these respective regimes the behavior of 
the entropy in the two theories is hardly comparable. 
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